Tuning the interactions of spin-polarized fermions using quasi-one-dimensional confinement 



m 
o 
o 



in 



> 

00 

in 
cn 

o 
m 
o 



o 



Brian E. Granger^ and D. Blume^ 

^ Institute for Theoretical Atomic, Molecular and Optical Physics, Harvard-Smithsonian CFA, Cambridge, MA 0213^ 
^Department of Physics, Washington State University, Pullman, WA 99164-281^ 

The behavior of ultracold atomic gases depends crucially on the two-body scattering properties of these 
systems. We develop a multichannel scattering theory for atom-atom collisions in quasi-one-dimensional (quasi- 
ID) geometries such as atomic waveguides or highly elongated traps. We apply our general framework to the 
low energy scattering of two spin-polarized fermions and show that tightly-confined fermions have infinitely 
strong interactions at a particular value of the 3D, free-space p-wave scattering volume. Moreover, we describe 
a mapping of this strongly interacting system of two quasi- ID fermions to a weakly interacting system of two 
ID bosons. 
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Two particle scattering is ubiquitous in physics. With the 
achievement of quantum degeneracy in ultracold atomic gases 
Ql^, renewed interest in the scattering of two atoms at low 
temperatures has arisen. This interest in atom-atom scattering 
stems from the fact that, to first order, the many-body physics 
of an ultracold atomic gas depends on a single atom-atom 
scattermg parameter For spin-polarized bosons, this pa- 
rameter is the 5-wave scattering length as — — limi^^o tan hs/k, 
where 5^ is the i-wave scattering phase shift at collision mo- 
mentum k. For two colliding spin-polarized fermions, in con- 
trast, the antisymmetric character of the wave function forbids 
i-wave scattering, and instead, p-wave scattering becomes 
dominant. The interaction is then characterized by the p-wave 
scattering volume Vp — —\\mk^Qt&nhp/lc' [4J. 

An exciting feature of ultracold atomic gases is that the ef- 
fective two-body scattering properties can be manipulated by 
external magnetic fieldsLSi or by strong confinement in 
one or more directions iTjIa]. Of particular interest are quasi- 
ID geometries created by atomic waveguides or optical lat- 
tices, which may, if loaded with ultracold bosons or fermions 
II9I I10II. pr ovide an opportunity to study novel ID many body 
states f 1 1], and to perform high precision measurements 1 12]. 
Because the many-body physics of such quasi- ID systems 
depends predominantly on the atom-atom scattering proper- 
ties in the waveguide, it is imperative to understand how the 
waveguide modifies the free space scattering properties. 

The scattering of two bosons in a quasi- ID geometry has 
been studied by Olshanii 1 13]. For a waveguide with harmonic 
confinement the 3D Hamiltonian for the relative coordinate 
r = (p,(|),z) of two atoms of mass m reads: 



(1) 



where fi^m/l denotes the reduced mass, and V3£)(r) the full 
atom-atom interaction potential. Using a regularized zero- 
range potential, Vioir) = ^h^asb{r)^r, Olshanii derives an 
effective ID Hamiltonian, 
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which reproduce the low energy scattering solutions of the 
full 3D Hamiltonian, Eq. Q. This elegant result ]13, 1% 
shows that the strong confinement of the waveguide gives rise 
to an effective ID interaction, parameterized by gw, which 
can be tuned to virtually any strength by changing the ratio 
of the 3D 5-wave scattering length a^ to the transverse os- 
cillator length a± — y/h/fj(0± through, for example, appli- 
cation of an external magnetic field. Notably, the effective 
ID interaction becomes infinitely strong (gio diverges) when 
the 3D scattering length takes on the particular ^«/fe value 
as/aj_ = 1/ /2)| « 0.6848, where ^(•) is the Riemann zeta 
function 11511 . 

The results obtained for quasi- ID bosons lead us to the 
following question: How do two spin-polarized /erm/oni un- 
der quasi- ID confinement behave? To answer this question 
we develop a general framework to obtain scattering solu- 
tions of the waveguide Hamiltonian, Eq. Q, that are appli- 
cable to bosons and fermions at any energy, as well as to a 
large class of quasi- ID scattering processes, single- or multi- 
channel in nature, at any energy. When applied to identical 
fermions, this framework shows that the confinement modi- 
fies the free-space scattering properties of two spin-polarized 
fermions significantly. In particular, we show that two spin- 
polarized fermions — just as two spin-polarized bosons — 
can have infinitely strong effective interactions for a finite 3D 
scattering volume Vp when confined to a quasi- ID geometry. 

A key quantity in our formalism is the familiar /T-matrix 
K^^{E) [161, which encapsulates the free-space (no con- 
finement) scattering physics of a given atom-atom potential 
Vj,D{r) at all energies. One advantage of using A'-matrices 
is that the complications of using zero-range potentials are 
avoided. For general multichannel collisions the 3D A'-matrix, 
K^^, can be determined accurately by a number of numerical 
techniques for a single channel scattering pro- 

cess it reduces to a diagonal matrix of scattering phase shifts 
tan5/. Below, we derive an effective ID ^T-matrix, K}^, that 
i) is written in terms of the known 3D ^T-matrix, K^^, and 
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ii) fully characterizes the effective interaction of two tightly- 
confined atoms. 

For now, we assume an interaction potential V3£)(r) such 
that the 3D K-matrix reaches its asymptotic limit at a distance 
Kc much smaller than the transverse confinement length ax 
{rc ^ ax)- While this condition may not necessarily be ful- 
filled for a realistic atom-atom potential, one can often con- 
struct a model potential that fulfills this requirement while 
mimicking the behavior of a realistic two-body potential as 
closely as possible. Configuration space can then be par- 
titioned into two regions: i) r < r^. The two-body poten- 
tial dominates, while the confining potential is negligible, ii) 
r > Tc- The confinement is felt, while the atom-atom potential 
is negligible. For r « r^., the full wavefunction can be conve- 
niently written in spherical coordinates 1 16|: 



(4) 



where the 3D /T-matrix, K^^, contains all the "scattering in- 
formation" about the atom-atom potential Vj,D{r). The en- 
ergy normalized regular solution f/„, (r) at energy E / h(0± = 
{kaj_)^/2 is given in terms of the spherical Bessel function 
jiikr): 



F,,„{r)=Yi„,{Q,<^) 
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The irregular solution Gii„{r) is identical to /v,„(r), but with 
ji{kr) replaced by ni{kr). In the spherical representation, the 
parity of the relative wavefunction is simply %tot — 

To obtain the ID scattering properties, the spherical solu- 
tion, Eq. (0}, must be propagated outward to large |z| where 
the cylindrical symmetry of the harmonic confinement dom- 
inates. For |z| ^ ax ^ r^-, the wavefunction is a product of 
2D harmonic oscillator wavefunctions <!>„„, (p,(|)) and free par- 
ticle solutions in z- We choose the energy normalized regular 
\(/„,„(r) and irregular X„,„(r) cylindrical solutions to be eigen- 
states of the parity 7t,„, = 7t,(— 1)'" : 



¥nm(r) 

Xnm(r) 



gn{z) 
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Here, the cylindrical channels are labeled by the harmonic os- 
cillator quantum number n and the angular momentum m, so 
that the z-direction momentum q„ in each channel is defined 
by the relationship. 



(fcax) 



ftcox 



= {2n + \ + \m\)- 



(7) 



The parity 7t- is determined by the forms of the regular /„ (z) 
and irregular gn{z) free particle wavefunctions. To always 
work with eigenstates of 71,. and to keep the regular and irreg- 
ular functions 90° out of phase at both positive and negative z, 
we use the conventions 12011 . 



cos{q„ \z\) for 71, = 1 
||Sin(^„|z|) for7i- = -1 ' 
(8) 



sm{q„ \z\ ) for7t, = 1 



Vtoxax V^^^ \ -^cos(^„|z|)for7i:, ==-1 



(9) 

The full asymptotic wavefunction at large |z| can then be writ- 
ten in terms of the solutions \(/„m(r) and Xn;n(r), and a ID 
K-msLtTix K^'^ S: 



n 



(10) 



To express the ID ^-matrix in terms of the 3D A'-matrix 
we use a frame transformation 121. 22,. 231 which transforms 
one set of solutions of the Schrodinger equation (labeled by a 
set of quantum numbers) to another set of solutions (labeled 
by a different set of quantum numbers). Such a transforma- 
tion is exact and orthogonal when the two sets of solutions 
satisfy the same Schrodinger equation everywhere. A local, 
nonorthogonal frame transformation, on the other hand, is 
useful for transforming between two sets of solutions that sat- 
isfy the same Schrodinger equation only in a limited region of 
space. Following Greene in his treatment of negative ion pho- 
todetachment in magnetic fields i24ll . we use the local frame 
transformation to relate our spherical free particle solutions, 
Eq. (0}, to our harmonically confined cylindrical solutions, 
Eq. (I10> . through a nonorthogonal matrix Uj. 

¥««i(r) = l^/7m(r)t//„ and x„,„{r) = ^G/„,(r) (t/^),^,' ■ 

(11) 

The sum over / in these expressions is understood to be over 
I — 0,2,4, . . . for 71, = 1 and over I — 1,3,5, . . . for 7t, — —1. 
The elements of the frame transformation matrix Ui„ are cal- 
culated by projecting the exp ressions given in Eq. il 1> onto 
the spherical harmonics [24]. Using the transformation ex- 
pressions, Eq. (EJj in Eq- < I10> . the ID /T-matrix K}'^ can be 
expressed in terms of K^^ at all energies: 



(12) 



This relationship is essentially exact for any Kj'^, including 
multichannel cases, as long as ^ a±. 

Up to now, the cylindrical asymptotic wavefunction ^a(r), 
Eq. (I10> . is written in terms of the regular and irregular cylin- 
drical functions solutions \|/nm(r) and Xnm(r), which contain 
exponentially diverging pieces in channels n that are energet- 
ically closed {E/h<X)± < 2n + 1 + \m\). To obtain asymptotic 
solutions with the correct, exponentially decaying boundary 
conditions in the closed channels, these divergences must be 
eliminated. To do this, we use the approach of multichannel 
quantum defect theory and partition K}^ into a closed ("c") 
and an open ("o") subspace 



part i 



KlS K 



(13) 



Eliminating the closed channels results in a "physical" K- 
matrix in the open channels, Kj,^''''"\ 



K 
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ID 
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ID 



(14) 
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The corresponding asymptotic wave function I'Jj'"' — a now 
denotes an open channel — having the correct physical bound- 
ary conditions in all channels involves a sum over only the 
open channels: 



<'"(r)= ^ <!>„ 

nGopen 



(PA) fn{z)?>na- gn{z)K, 



lD,phys 
na 



(15) 

Equations (I12> . (I14> and il5\ provide a rigorous path from the 
full 3D scattering properties encapsulated in K^^{E), to an ef- 
fective ID system, Eq. M5\ . whose scattering properties are 
given by K^Jo'''^^\E). Our framework shows that the closed 
channels contribute significantly to the effective ID scatter- 
ing properties at energies E near closed channel resonances, 
where det[i - iK^^ (£)] w (see alse Q). 

We now introduce an important simplification. When the 
3D scattering properties are dominated by the phase shift 
tan 5/ in a single partial wave / (Kp^, = 5;";' 5//' tan 5/), the ID 
/T-matrix, Eq. (I12> . becomes a rank one matrix, 

,tan biUi„t. 



K. 



ID 



(16) 



Thus, the closed channel part of K}^ also has rank one with a 
single eigenvalue Xc'. 



X,=Tr^.f =tan5, ^ (f/„,)' 

riGclosed 



(17) 



This fact allows Koo'''''^'\ Eq. il4l . to be simplified by diago- 
nalizing and inverting the matrix [i — analytically: 

(18) 

Equation ( I18t shows that the "bare" ID A'-matrix K}^^ is 
renormalized by the closed channels physics encapsulated in 
the eigenvalue Xc{E), which can be strongly energy depen- 
dent. As the relative energy E increases, new channels be- 
come open and the dimension of Koo''^'^-'^ increases to reflect 
the multichannel nature of the scattering. 

Application of our framework to the scattering of two iden- 
tical bosons under strong transverse confinement predicts the 
following general result: when as{E)/a± = 0.6848, a diver- 
gence of the effective ID interaction occurs at all threshold 
energies, E = h(x)± (2n + 1 ) for m = 0, where a new transverse 
mode becomes open. This extension of Olshanii's zero energy 
result 1 13| to all energies and multichannel cases suggests that 
a Tonks-like regime |J_3] can be realized experimentally in a 
"high-temperature" gas, ksT 3> ftcox- Further details regard- 
ing our results for bosons will be published elsewhere '25^. 

We now present an application of our framework to the 
low energy scattering (m = 0, 1 < E/h(0± < 3) of two spin- 
polarized fermions, whose 3D scattering properties are param- 
eterized by a single parameter, the energy dependent scatter- 
ing volume Vp {E) = — tan 5,, (£) /k^ . In this case, the resulting 
effective fermionic ID /^-matrix, Eq. J18t . becomes: 



K 



w,- 



K 



W.phys _ 



6V„ 



-qoa± 



(19) 



Here, the Hurwitz zeta function ^(-j) lIlSll arises from the 
eigenvalue Xc of K}.^, Xc — Tr^.^. The /T-matrix given in 
Eq. ( I19I I along with the odd-parity wavefunction. 



¥ W^uT sin(^ok[ 



cos 



10 \z\ )Kqq' 



(20) 



provides a complete scattering solution to the waveguide 
Hamiltonian, Eq. Q, when a single cylindrical channel is en- 
ergetically open. The ID /T-matrix, Eq. il9i . diverges when 
the scattering volume Vp has the particular value 



ycnt 



12C 



1 1 

"2'2 



(21) 



This implies that two spin-polarized quasi- ID fermions have 
infinitely strong interactions for a finite 3D scattering volume 



Vn 



ycnt 
"p ■ 



Next, we derive an effective ID Hamiltonian that describes 
many of the low energy properties of two spin-polarized 
fermions in a waveguide [see Eq. ^]. Importantly, the ID 
zero-range potential giD?>{z) [see Ecl (01, which has been 
very successful in treating bosons II ll 11311 . cannot be used 
directly since it results in an unphysical scattering amplitude 
for fermions. One way around this difficulty would be to 
use a zero-range potential that gives a meaningful scattering 
amplitude for fermions |26|. Alternatively, we propose to 
map the fermionic K^^-^ , Eq. M9\ . to a bosonic ID K-matrix 
(along with the corresponding wavefunctions). Mappings be- 
tween fermions and bosons are important in theoretical treat- 
ments of ID many body systems, as they allow one to under- 
stand systems of strongly interacting ID bosons (fermions) by 
mapping them to weakly interacting systems of ID fermions 
(bosons) U_l]. 

At low energies, the ID scattering wavefunction \|/^ {£) for 
two fermions is given by Eq. j20> while that for two bosons 
reads: 



\(/+(z)- [cos(^ok|)-sin(^ok|)^^ 



io,+l 



(22) 



where + denotes the even parity ID /T-matrix. With the 
choice. 



K, 



\D. 



ID. 



(23) 



^00 ~ ^/■'^oo ' 

the bosonic wavefunction can be written in terms of the 
fermionic one (and vice versa): 



(24) 



Application of the proposed mapping, Eqs. \23\ and ( I24> . to 
our effective ID /T-matrix for two fermions, Eq. (I19> . results 
in an equivalent system of two ID bosons interacting through 
the potential §jq''5(z), with the "mapped coupling constant": 



map 
8lD 

flxfi-K>I 



,3 r 
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a^ 



1 1 
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(25) 
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This remarkable result implies that two spin-polarized quasi- 
ID fermions with infinitely strong interactions [see Eq. ( I21H . 
can be mapped to a system of non-interacting bosons. More 
importantly, however, the mapped g'"^'', Eq. i25\ . applies to 
any interaction strength Vp. 

To confirm our analytical results we have performed numer- 
ical calculations for two spin-polarized fermions in a highly 
elongated trap interacting through a two-body model poten- 
tial, Vioir) = d / cosh^ {r / b) . We use a B-spline basis set 
to solve the 3D Schrodinger equation for the Hamiltonian 
given in Eq. Q with the confining potential replaced by 
5/jC0^(p^ + A,^z^), where A, = 0.01. The antisymmetry of 
the fermionic wave functions is enforced by only including 
odd-parity states. Solid lines in Fig.^show the resulting en- 
ergy spectrum, E > h(0±, as a function of the two-body well 
depth d. Straight dotted lines indicate the energy levels of two 
non-interacting spin-polarized fermions, while straight dashed 
lines indicate those for two non-interacting bosons. As pre- 
dicted analytically (and indicated by an arrow in Fig. [0, the 
full 3D energy levels coincide with that of two non-interacting 
bosons for well depths d corresponding to Vp = Vp"' . 

FigurefOadditionally shows the spectrum of the ID bosonic 
Hamiltonian, Eq. (|2}, with the additional potential i/jA,^co^z^, 
using the mapped coupling constant g'"^'' [Eq. (I25H (aster- 
isks). To account for the energy-dependence of §'"d'' 
consistently solve for the scatterin g vo lume Vp and the ID 
energy levels for each well depth d 12 711 : details will be pub- 
lished elsewhere [25]. Figure [0 shows that the ID energies 
agree very well with the 3D energies. Two interacting spin- 
polarized fermions under quasi- ID confinement can hence be 
mapped to a system of two ID bosons interacting through a 
5-function potential with a mapped coupling strength g"'^'', 
Eq. ( I25> . We note that present-day experiments |6] can ac- 
cess regimes where our theory applies; the question of how to 
experimentally verify our predictions, however, is beyond the 
scope of this paper. 

In summary, this Letter develops a general framework for 
treating atom-atom scattering under quasi- ID confinement. 
The framework includes multichannel collisions at arbitrary 
energies and is limited only in the assumption that the char- 
acteristic length of the two-body potential is small compared 
to the confinement length a±. Application of the framework 
to two spin-polarized fermions in a waveguide shows that the 
3D scattering behavior, which is assumed to be characterized 
completely by the scattering volume Vp, is altered signifi- 
cantly due to the presence of the tight transverse confinement. 
Specifically, two quasi- ID fermions can have infinitely strong 
interactions for a finite 3D scattering volume. Finally, the sys- 
tem of two quasi- ID fermions can be mapped to a system of 
two ID bosons. 

Further study is required to understand the full implica- 
tions of these results for the many body physics of quasi-lD 
spin-polarized fermions. At this point, we only mention that 
any many body theory will have to carefully account for the 
energy-dependence of the mapped coupling constant ^i'd''- 

Discussions with Chris Greene and Mike Moore were help- 
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Figure 1: 3D energy levels (solid lines) together with ID energy lev- 
els (asterisks) as a function of the two-body well depth d (b is fixed at 
b = 0. 13ax) for two spin-polarized fermions under highly-elongated 
confinement. Straight dotted lines indicate the energy levels of two 
non-interacting spin-polarized fermions, while straight dashed lines 
indicate those for two non-interacting bosons. Arrows indicate val- 
ues of d for which the corresponding scattering volume Vp (E) has 
the particular values Vp = and Vp = V"" , respectively. 
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